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Abstract 


The object of the paper is to develop a simple and practical control strategy for a molten carbonate fuel cell (MCFC) system. Based on a dynamic 
model and the control demand, a cascade control strategy is designed. The master controller imposes a change of cell current representing a change 
in power demand and sets the amount of fuel gas, the steam-to-carbon ratio, the air number and the cathode gas recycle ratio to their corresponding 
conditions for optimal steady state electric efficiency. Two feedback PID controllers are in the inner loop, one guarantees the solid temperature not 
to exceed a maximum temperature by changing the air number around the default set by the master controller; the other controls the maximum 
temperature difference by adjusting the steam-to-carbon ratio. Step response tests show that this control strategy works well when the cell current 


changes. 
© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


Molten carbonate fuel cells (MCFC) are a moderate to 
high temperature fuel cell technology. The nominal operating 
temperature is around 600°C, which is the typical operating 
temperature found in conventional utility power plants. The- 
oretically, the MCFC has a high fuel-to-electricity efficiency 
(~55-60%) and minimal environmental emissions such as nitro- 
gen oxides (NOx) (<1 ppm). These properties make the MCFC 
an attractive candidate for decentralized power generation 
[1]. 

An industrial MCFC stack is characterized by a high degree 
of system integration, as the process behavior depends on 
numerous interactions between the electrochemical reaction 
steps, the internal reforming, mass transport processes, and the 
heat transfer inside a cell or stack. An understanding of the 
physical—chemical processes in an MCFC can be obtained from 
mathematical models based on physical conservation principles 
[2]. Model-based process control and process design strategies 
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can lead to a much better use of the fuel cells’ capacities and 
increase the efficiency of the system. 

As the overall fuel-to-electricity conversion efficiency is the 
most concerned interest for power industries, a maximum elec- 
tric efficiency is desired not only when the cell is operated at 
constant load over a long period of time but also during dynamic 
load changes. In the former situation, an optimal electric effi- 
ciency can be achieved through a steady state optimization pro- 
cedure [2]; whilst in the latter situation a fast and safe transition 
must also be taken into consideration. To keep the process safe, 
the solid temperature must not exceed a minimum and a maxi- 
mum temperature, and also the maximum temperature difference 
within the solid phase is limited. In literature, only a few publica- 
tions about model-based control of MCFC are available. Lucas 
and Lee [3] describe a first principle-based non-linear dynamic 
model and also obtain reduced-order models, but no spatial gra- 
dients are considered. Kang et al. [4] derive a linear 3 x 3 transfer 
function matrix model from experimental study, which will sim- 
plify the design of controllers but will limit the control is to the 
vicinity of a chosen operating point. Shen et al. [5] present an 
adaptive fuzzy control procedure for the temperature of MCFC 
stack based on a neural network identification model. As the 
dynamic model is also empirical, the applicability of the con- 
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Nomenclature 

Teel total cell current 

S/C steam-to-carbon ratio 

S/C steam-to-carbon ratio estimated from steady state 
optimization 


Ucell cell voltage 


Greek symbols 

Ifeea amount of fuel gas 

Nel electric efficiency optimized 
Hel electric efficiency controlled 


D gift maximum temperature difference 

Omax maximum temperature 

Dmin minimum temperature 

Xair air number and cathode gas recycle ratio 

air air number and cathode gas recycle ratio 
estimated from steady state optimization 


troller is restricted to operation conditions, where this model is 
valid. This work tries to construct a framework for satisfying 
high electric efficiency over a larger range of operation points 
and at the same time controlling the solid temperature when 
power load changes. For simulation, optimization and control of 
MCEC processes mathematical models are required. Firstly, the 
paper summarizes a detailed spatially two-dimensional model of 
the MCFC presented in previous publications [2,6]. A control 
strategy is then constructed which achieves the optimal electric 
efficiency by a steady state optimization and at the same time 
satisfies the solid temperature constraints by involving two PID 
controllers. 

Currently, the combination of independently designed feed- 
forward and feedback controllers [12] is attracting increasing 
attention as a control concept for nonlinear systems [13]. Aguiar 
etal. [14] use such a concept for the control of the outlet fuel tem- 
perature of a SOFC and test the controller in simulations with a 
spatially one-dimensional model. Spatial temperature gradients 
are monitored in [14], but it is not tried to control them directly. 

As will be shown in the later sections, the control approach 
proposed in this paper includes the control of spatial dependen- 
cies (the control of stack temperature difference). It combines 
a feed-forward branch based on a detailed physical model and 
feedback PID controllers. In practice, a suitable state estimator 
as presented in [6,10] will make the control method applicable, 
if the maximum temperature and the maximum spatial temper- 
ature difference cannot be measured. 

The application example considered in this paper is the so- 
called HotModule by MTU CFC Solutions [7]. The HotModule 
is an integrated MCFC system with a cell stack of 342 fuel cells 
that can deliver up to 250 kW of electrical power. A pilot plant 
installation of the HotModule exists at the Magdeburg university 
hospital, Germany. This pilot plant has been subject of a joint 
research project of several academic and industrial partners that 
aimed at model development and model-based process control 
[6]. Within this project, the HotModule has been used for experi- 


mental validation of the developed physical models [2] and state 
estimators [6,10]. The controller presented here is an extension 
of the results obtained in the project. 


2. Dynamic model 
2.1. Working principle of an MCFC 


Fig. 1(a) shows the principle of an MCFC with direct internal 
reforming (DIR). The MCFC consists of two porous electrodes, 
an electrolyte layer between them and gas channels above each 
electrode. The anode channel is fed with a mixture of natural gas 
(methane) and steam. In an internal reforming step, the reactions 
are oxidized to hydrogen and carbon monoxide. Both are oxi- 
dized electrochemically at the anode, consuming carbonate ions 
from the electrolyte and producing free electrons. 

On the cathode side, new carbonate ions are formed in a 
reduction reaction under consumption of electrons from the elec- 
trode. The electrons from the anode are transferred to the cathode 
via an external electric circuit, where they perform electric work. 

In this work, a MCFC system is considered, whose simplified 
scheme is depicted in Fig. 1(b). The fresh anode feed consisting 
of steam and methane first enters a pre-reformer where a part of 
the methane is converted to hydrogen, and then is distributed to 
the anode gas channels of the cell stack. The outlet flow on anode 
side goes into a burner. The burner is supplied with air from the 
ambient and oxidizes unconsumed methane and hydrogen. The 
burner outlet is fed to the cathode gas channels. Finally, a part 
of the cathodic exhaust gases may be recycled to the burner. 


2.2. Two-dimensional spatially distributed MCFC model 


A detailed model of a planar fuel cell with cross flow has 
been derived from first principles [8,9]. The model is dynamic 
and spatially distributed in two space coordinates. The main 
model assumptions are: 


e Isobaric gas phases, ideal gas law. 

e Plug flow in all channels. 

e Allimmobile parts of the cell, i.e. electrodes, electrolyte, bipo- 
lar plates, etc. are lumped into a solid phase with respect to 
their enthalpy balance. 

e Reversible potential kinetics for reforming reactions, 
Butler—Volmer kinetics for electrochemical reactions. All 
kinetics feature an Arrhenius term. 

e The electric potential field is spatially distributed and allows 
for the calculation of an inhomogeneous current density dis- 
tribution. 

e Mass transport resistance in the electrode pores is considered. 


The resulting model is a system of hyperbolic and parabolic 
partial differential equations for the temperatures in the gas 
and solid phases, for the composition of the gases, and for the 
electrical potentials. A reduced version of the model, which is 
more suitable for process control purposes, has been derived in 
[10,11]. 
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Fig. 1. (a) Working principle of an MCFC with direct internal reforming (DIR) and (b) scheme of the MCFC system considered in this work, which consists of a 


pre-reformer, the cell stack, and a burner between anode side and cathode side. 


3. The control strategy 


The model described in Section 2 gives the basis to design a 
control strategy. A fast and safe dynamic response during the cell 
current change is necessary. As the power load does not change 
very frequently an optimal steady state electric efficiency is also 
desired. In the ideal case the electric efficiency is dynamically 
optimized such that both demands can be achieved at the same 
time. However, the optimization of a complex non-linear fuel cell 
model in real time is very difficult. To overcome this problem, 
the approach adopted here is to achieve the optimal state electric 
efficiency in offline optimization and then to make sure that the 
dynamic response is as fast and safe as possible. That is, when 
the power load changes (i.e. the cell current changes) a master 
loop switches the operating variables to their new static condi- 
tions that satisfy an optimal electric efficiency at the new load. 
Then in the inner loop feedback controllers satisfy the dynamic 
performance criteria. In the following, steady state optimization 
will be introduced first and temperature response tests will then 
be performed in order to decide which operating variable con- 
trols which controlled variable. Control schematic and results 
will be given thereafter. 


This approach is similar to the control concept proposed 
by Aguiar et al. for an SOFC stack [14]. However, in con- 
trast to [14], a spatially two-dimensional temperature field is 
considered in this work. Furthermore, the concept proposed 
here controls the maximum temperature and the maximum spa- 
tial temperature difference at the same time. This leads to a 
multiple-input/multiple-output (MIMO) control problem, where 
the choice of suitable manipulated variables is more difficult than 
in the single-input/single-output case. 


3.1. Steady state optimization 


Heidebrecht [2] optimizes the input parameters of a single 
cross flow cell at given cell current to yield an optimal electric 
efficiency while fulfilling several temperature restrictions and 
avoiding carbonization in the anode channel. The optimization 
is performed at different cell currents, so that several points of 
a current-voltage curve with optimal operating conditions are 
evaluated. 

The electric efficiency, which is used as the objective func- 
tion of the optimization relates the total electric cell power 
diminished by the system’s parasitic power consumption to the 
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combustion enthalpy input of the fuel gas. The optimization 
has several constraints concerning the solid temperature and the 
possibility of carbonization. The temperature must not exceed 
a minimum and a maximum temperature, and also the maxi- 
mum temperature difference within the solid phase is limited. 
Three operation parameters, which can be adjusted easily in the 
HotModule, are chosen as optimization variables (cf. Fig. 1b). 
The first optimization variable is the total amount of fuel gas 
(teed) fed to the system. In the case of the HotModule, the 
fuel gas consists of steam and methane. The ratio between the 
flux of steam and the flux of methane fed to the system, i.e. 
the steam-to-carbon ratio (S/C), is chosen as the second opti- 
mization variable. It should be noted that water is generated 
in the anodic reaction and consumed in the internal reform- 
ing step. Therefore, the amount of water in the system affects 
the reforming reaction as well as the electrochemical reaction. 
The interaction between the exothermic anodic reaction and the 
endothermic reforming step strongly influences the temperature 
in the cell. Therefore, the steam-to-carbon ratio has a strong 
impact on the overall system behavior and can be seen as one of 
the most important operation parameters. This is also reflected 
by the fact that in the case of the HotModule the steam-to-carbon 
ratio can be set directly at the process control system. The third 
optimization variable is the air number (A,ir), defined as the ratio 
between the amount of air actually fed to the burner and the 
amount of air required under stoichiometric conditions. Those 
three variables must be adjusted to meet the power demands 
from an external load on the cell in an optimal way. For sim- 
plicity, the effect of the external load is modeled by setting the 
total cell current (Ice) to a fixed value within each optimiza- 
tion. 

Some typical input quantities at the optimum points together 
with additional information like temperature are listed in Table 1. 
For the control purpose the number of steady state optimization 
points should be infinite to obtain optimal steady state perfor- 
mance at arbitrary cell currents. In order to have more possible 
operation points the table by Heidebrecht [2] was extended by 
specifying Icey change by 0.01 from 0.5 to 1.0. When an Jey 
change is less than 0.01, the optimized variable values between 
two nearest points are interpolated. For an Kcen change, the mas- 
ter controller looks up the table, locates the column of the desired 
Icey, State and sets the three variables (I feeg, S/C, Aair) to their 
corresponding numbers. 


Table 1 

Results of the steady state optimization of input parameters 

Teei 0.5 0.6 0.7 0.8 0.9 1.0 

T feed 0.6783 0.8296 0.9841 1.114 1.2859 1.2920 
S/C 2.4304 2.4707 2.4903 2.4321 2.3106 1.7049 
air 2.0250 2.0659 2.1913 2.3447 2.6071 2.9574 
Din 2.9604 2.9999 2.9999 2.9919 2.9942 3.0000 
O max 3.1604 3.2000 3.1999 3.1994 3.2017 3.2000 
Date 0.2000 0.2001 0.2001 0.2075 0.2076 0.2000 
Nel 0.5391 0.5240 0.5030 0.4809 0.4384 0.3838 
Mel 0.5391 0.5240 0.5030 0.4769 0.4352 0.3838 


jel refers to the closed-loop case discussed in Section 3.4. 


3.2. Solid temperature response analysis 


The control object is to fulfill the temperature constraints 
when the current load changes, that is, to transfer the fuel cell 
safely from an optimal steady state to another. The steady state 
optimization introduced in the previous section gives a reference 
of what the input parameters (I feea, S/C and Agir) would be like 
to achieve the desired new state, but the dynamics of transition 
do not guarantee the temperature constraints. In order to control 
the solid phase temperature during dynamic transition, tempera- 
ture control loops are required. The possible controlled variables 
are minimum temperature (J min), maximum temperature (max) 
and maximum temperature difference (vit). It is found that the 
minimum temperature rarely exceeds the limits, so there are two 
controlled variables (max and Ŷqiff) and as stated already three 
manipulated variables (I"feea, S/C and Aair). This is a typical 
multi-input—multi-output (MIMO) control problem. In the fol- 
lowing, a match of input and output variables based on physical 
considerations is attempted. Dynamic responses of the system 
model to step changes of the three manipulated variables are 
considered. Fig. 2 shows the step responses of all the possible 
input and output matches when each of the three input variables 
changes from its corresponding optimized state at Icey) =0.8 to 
the optimized state at Ke1 =0.7 at time 5000 and then changes 
back to the optimized state at Jeen =0.8 at time 10,000. The 
actual cell current is kept constant at ee =0.8. All results are 
given in dimensionless numbers, one dimensionless time unit 
corresponding to 66s. For the cell current and the temperature, 
the following scaling is used: one current density unit equal to 
840 A m7”, one temperature unit equal to 273 K. 

Fig. 2(a) gives the temperature difference and maximum tem- 
perature response when switching Ifeeq from its optimal state 
(1.114) to the new optimal state (0.9841) while keeping Air and 
S/C at their old state. A reduction of I"feeq means physically that 
less fuel is fed to the system. As the cell current is kept constant 
and the cell voltage changes only slightly, the electrical power 
hardly varies during the simulation experiment in Fig. 2(a). Con- 
sequently, the reduction of Ifeeq results in a higher electrical 
efficiency and — because less chemical energy is converted to 
heat — in a lower average temperature of the stack. However, 
the new steady state with the higher efficiency, which is reached 
at time ~6000, is not feasible as the temperature difference is 
above the tolerable limit. This behavior is found to be typical 
for changes of Ifeeq at different operation points. A reduction 
of Ife by a small amount always leads to a spatially more 
inhomogeneous temperature profile and hence to an increase of 
the maximum spatial temperature difference. The implications 
of a modified Ifeeq on the absolute maximum temperature are 
less clear. Depending on the actual operation point, the temper- 
ature maximum may increase or decrease. Therefore, Ifeeq iS a 
suitable variable for controlling the temperature difference but 
is not suitable for controlling the maximum temperature. 

The system response to a variation of the steam-to-carbon 
ratio is qualitatively similar to the response to a variation of 
Iteea. An Example is shown in Fig. 2(b). An increase of S/C 
at a constant J"feeq also means a reduction of the amount of 
fuel fed to the system. Therefore, the argumentation used to 
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Fig. 2. (a) Omax, Vaife response to a I feeq change; Ifeeq changes at from the optimal value at Ke1 =0.8 to the optimal value at Ke1 =0.7 and back again. (b) Umax, 
aif response to a S/C change; S/C changes at from the optimal value at Fee =0.8 to the optimal value at Zei =0.7 and back again. (c) max, aiff response to a Aair 
change; Air changes at from the optimal value at ce, = 0.8 to the optimal value at J¢e = 0.7 and back again. 


explain Fig. 2(a) holds for (b), too. An increase of S/C generally 
causes an increase of the temperature difference, while the max- 
imum temperature may decrease or increase, depending on the 
operation point. Therefore, the steam-to-carbon-ratio is another 
possible candidate for the control of the spatial temperature dif- 
ference. 

The effect of a varying air number air is shown in Fig. 2(c). A 
higher amount of air in the system mainly has a cooling effect, i.e. 
it reduces the maximum temperature. The influence on the spa- 
tial temperature difference is not that strong. Therefore, varying 
Aair locally around a given set-point is a good method to control 
the maximum temperature of the fuel cell. 

To summarize the discussion of Fig. 2, the air number Aair has 
the most direct relationship to max because of its cooling effect. 
As to Vaite, the effect of Ifeeq and S/C is almost the same and in 
this work S/C is chosen to control Maire. The resulting cascade 


control scheme is depicted in Fig. 3. A feed-forward controller 
sets the three manipulated variables Ifeea, air, S/C to the opti- 
mal steady state values Ifeeq, Ñair, S/C. The air number Air 
is tuned by the PID controller 1 around the default set by the 
feed-forward controller in order to obtain a pre-defined max- 
imum solid temperature max also under transient conditions; 
similarly, S/C is tuned by PID controller 2 around the optimum 
steady state value in order to achieve a pre-defined maximum 
temperature difference V gir. 


3.3. Control results and analysis 


Fig. 4(a) and (b) present the ?gig¢ and Pmax transient behavior, 
under open-loop and closed-loop conditions, when the cell cur- 
rent changes from 0.55 to 0.6 and further to 0.65, 0.75, 0.8, 0.85, 
0.9, 0.95, 1.0 at every 2000 time units (0, 2000, . . ., 18,000) rep- 
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Fig. 3. The cascade control scheme. 


resenting positive load changes and then changes back from 1.0 
to 0.95, 0.9, ..., 0.5 at times (20,000, 22,000, . . ., 38,000) rep- 
resenting negative load changes. Note that closed-loop implies 
that all control loops are active, whereas open-loop suggests that 
only the master controller is in operation. 

The following parameters of the PID controllers are used: 
Vaig controller gain=180.0, integral time=20.0, derivative 
time = 500.0; max controller gain = 24.0, integral time = 0.29, 
derivative time = 40.0. These parameter values are found to pro- 
vide a satisfactory system response. As can be seen, for the 
open-loop case (no temperature control), Vgir¢ and max usually 
exceed their corresponding limits (0.2 and 3.2, respectively), 
while for closed-loop case, the PID controllers successfully and 
quickly bring the temperatures back to their desired set-points. 


T_diff 


T_MAX 


(b) TIME x 10° 
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In the simulation shown in Fig. 4, the desired aire was always 
set to 0.2, resulting in a steady state electric efficiency deviating 
slightly from the optimal value, as can be seen from Table 1. 


3.4. State estimator 


Usually, the required measurement information on v gifs and 
Vmax 1s not directly available from an industrial fuel cell system. 
In most cases, temperature measurements can be taken only at 
a few points in the cell stack. To make the proposed control 
scheme applicable in practice, a state estimator or observer is 
required. Such state estimation techniques were developed in 
previous work [6,10] based on a reduced MCFC model [11]. It 
could be shown that the developed observers provide accurate 


a) 1 2 3 4 
TIME x10" 


Fig. 4. (a) aire control results (open-loop on left panel and closed-loop on right panel) and (b) max control results (open-loop on left panel and closed-loop on right 


panel). 
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estimates of the spatial temperature profile, but use only a small 
amount of CPU time. 


4. Conclusions and outlook 


In open literature, comparatively few publications are avail- 
able on model-based process control of MCFC systems. The 
results of this work show that a rather simple control scheme can 
contribute substantially to a safe and efficient operation of a high 
temperature fuel cell stack. Such systems are rather sensitive 
to local over-temperatures and spatial temperature differences. 
Therefore, temperature management is a key problem of the pro- 
cess operation of high temperature fuel cells. This work proposes 
a cascade control scheme that may overcome this problem. The 
cascade control consists of three control loops. The first loop 
is composed of a master (feed-forward) controller that imposes 
a load change and sets the fuel gas, the steam-to-carbon ratio, 
and the air number to their corresponding conditions for optimal 
steady state electric efficiency. The master-controller reads the 
steady state operation conditions from a look-up table. In this 
work, the look-up table was generated by optimizing the steady 
state solutions of a detailed process model; alternatively, it would 
also be possible to determine optimal steady state operation con- 
ditions experimentally and store them in the look-up table. 

The other two loops are feedback PID controllers that are 
responsible for the temperature management under transient 
conditions. One of the PID controllers adjusts the air number 
in order to keep the maximum temperature below a pre-defined 
limit. The other PID controller corrects the steam-to-carbon ratio 
in order to limit the maximum spatial temperature difference. 
This simple control structure with linear feed-back controllers 
is inspired by the physical properties of the system and offers 
several advantages. Firstly, the actions of the controller appear 
to be transparent and understandable to the operating personal, 
a demand often made by practitioners on advanced control 
schemes. Secondly, the scheme has been derived from rather 
general properties of high temperature fuel cells and, therefore, 
may be robust to model uncertainties and applicable to many 
fuel cell systems. 

The implementation of the developed cascade control is espe- 
cially simple, if it is possible to measure directly the maximum 
temperature and the maximum temperature difference, as was 
assumed in this work. In practice, this measurement information 
is not always available. Very often, temperature measurements 


can only be taken at the outlets of the gas channels or at a few 
distinct points within the cell stack. In such cases, it possible to 
estimate the desired maximum temperatures and temperature 
differences with good accuracy from the available measure- 
ment information by using state estimation techniques. This was 
demonstrated in previous publications [6,10]. In a next step, the 
controller developed here will be combined with a suitable state 
estimator in order to make the scheme applicable under condi- 
tions, where measurement information is limited. 
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